The stability and free expansion of a dipolar Fermi gas 
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We investigate the stability and the free expansion of a trapped dipolar Fermi gas. We show that 
stabilizing the system relying on tuning the trap geometry is generally inefficient. We further show 
that the expanded density profile always gets stretched along the attractive direction of dipolar 
interaction. We also point out that by switching off the dipolar interaction simultaneously with the 
trapping potential, the deformation of momentum distribution can be directly observed. 
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Introduction. — Since the predication of dipolar ef- 
fects in a Bose-Einstein condensate (l|, dipolar conden- 
sates have drawn significant interests during the past few 
years [2] . The anisotropic nature and the associated tun- 
ability of dipole-dipole interaction provide us an unique 
platform to study the novel quantum phenomena. In par- 
ticular, shortly after the experimental realization of the 
dipolar condensate of Cr atoms |3|] , experimentalists now 
have control over both the contact interaction and the 
trap geometry, and have demonstrated the collapse and 
stabilization of the Cr condensate [3] ■ 

The theoretical study on ultra-cold dipolar Fermi gases 
was initiated by Goral et al. [5|. They analyzed the 
equilibrium state properties and its stability based on 
a variational approach. The subsequential work includes 
investigating the BCS pairing originating from the at- 
tractive dipolar interaction 'G\ and the strongly corre- 
lated states of rotating dipolar Fermi gases |7J] . In gen- 
eral, for fermionic atoms, the Pauli exclusion may com- 
pletely wash out the dipolar effects originating from the 
weak magnetic dipole-dipole interaction [5[. However, 
with the fast development in making ultra-cold polar 
molecules Q, we expect the dipole-dipole interaction to 
play an important role in determining the fundamental 
properties of the ultra-cold fermionic polar molecules. 

In this paper, we investigate the equilibrium state 
properties and the free expansion of a polarized dipo- 
lar Fermi gas. Based on the full numerical calculations, 
we show that the critical dipolar interaction strength is 
significantly lower than that predicted variationally, and 
it has a very weak dependence on the geometry of the 
trapping potential. We therefore propose to stabilize 
a dipolar Fermi gas by tuning the dipolar interaction 
strength via a rotating orienting field. We then show 
that, independent of trap geometry, the Fermi gas is al- 
ways stretched along the attractive direction of dipolar 
interaction during free expansion. Similar to the dipolar 
condensates [3|, |9|, [lOJ , the anisotropic expansion can be 
used as a diagnostic tool for dipole-dipole interaction. 

During the preparation of the manuscript, we become 
aware of a preprint by Miyakawa et al. [Ill ] . They revis- 
ited the equilibrium state properties of a dipolar Fermi 
gas using a variational ansatz capable of describing the 
deformation in momentum space. It was shown that the 



exchange dipolar interaction induced the deformation in 
momentum distribution and lowered the stability of the 
dipolar Fermi gas. We remark that in the treatment pre- 
sented in this work, the momentum distribution is as- 
sumed to be isotropic, and consequently, the exchange 
interaction plays no role in determining the equilibrium 
density and the dynamics. We shall determine the valid- 
ity region of our work and point out that the deformation 
in momentum distribution can be measured experimen- 
tally by switching off the dipolar interaction during free 
expansion. 

Model. — We consider a system of N fermionic po- 
lar molecules with permanent dipole moment d at zero 
temperature. For simplicity, we assume that all dipoles 
are polarized by an external electric field E which forms 
an angle tp to z-axis. To make the dipolar interaction 
tunable, we further assume that the orienting electric 
field fast rotates around z-axis. The time average of the 
dipole-dipole interaction potential becomes [13 | 
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where Cd = rjd'^ /{Atteo) with i] — {3 cos2 (p — l)/2 being a 
parameter continuously tunable within the range [—21 1] • 
The inclusion of 77 not only makes the sign of dipolar 
interaction changeable, it can also be used to completely 
switch off the dipolar interaction if ip equals to 54.7°, the 
'magic angle'. 

Within the semi-classical framework, the state of the 
system is described by the phase space distribution func- 
tion f{r,v,t) which, in the collisionless regime, satisfies 
the Boltzman-Vlasov kinetic equation [lj| 
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is the trapping po- 



tential and for simplicity, we assume that C/ho is ax- 
ially symmetric with tOx = Wy — i^±, C/dd(r) = 
Cd J dr'Vi{r — r')ri(r') is the Hartree-Fock mean-field 
term originating from the dipole-dipole interaction, and 
n{r) — J dvf^TjVjt) is the density distribution function. 
Since directly solving Eq. Q is generally inapplicable. 
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FIG. 1: The dipolar interaction dependence of the density 
deformation (a) and the volume of the cloud (b) for A = 0.1 
(dotted line), 0.7 (solid line), 1 (dashed line), and 10 (dash- 
dotted line). 



we shall make use of the scaling ansatz defined as 

/(r,v,t) = /o(R(i),V(i)), 



(3) 



where Rj{t) — rj/bj{t), Vj{t) — bj{t)vj — bj{t)rj, and 
/o is the equilibrium phase space distribution function. 
This scaling transformation was first introduced to s tudy 
the free expansion of a Bose-Einstein condensate 14 1, 
and the collective oscillations of a classical Bose gas 15 1 . 
Recently, it was also generalized to study the free ex- 
pansion and the collective excitation of both normal and 
superfluid Fermi gases [ij, [3 . Following the standard 
procedure [l3,[iMllg|j the Boltzman-Vlasov equation ^ 
can be reduced to the coupled equations for the scaling 
parameters bj, which, in dimcnsionlcss form 17|, reads 
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where \j = ljj/lj^, D = i^N^/^^/rnFuTi/W^^ is a di- 
mensionless quantity characterizing the strength of dipo- 
lar interaction, and (i?^) = J dIlR^no(R.) with no(R) = 
/ dV/o (R, V) being the equilibrium density. For the 
axially symmetric trap studied in this work, we have 
(Rl) = {R^)- Furthermore, 
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FIG. 2: The A dependence of the critical dipolar interaction 
strength D\ > (solid line) and D*_ < (dashed line). 



is the dipole-dipole interaction potential under the scal- 
ing transformation. Once the equilibrium density no is 
obtained, the dynamics of the system can be studied by 
evolving Eqs. ([4]). Before presenting our results on free 
expansion, we would like to address the equilibrium state 
properties of a trapped dipolar Fermi gas. 

Equilibrium density. — Under local density approxi- 
mation, the equilibrium density is determined by equa- 
tion [| 
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where A = A^ and the chemical potential /i is introduced 
to ensure that no is normalized to unit. The properties of 
equilibrium density was previously analyzed by Goral et 
al. [5| using variational approaches. Here we would like 
to tackle the problem numerically. Utilizing the cylindri- 
cal symmetry of the system and the numerical technique 
developed by Ronen et al. [l8|, we are able to obtain no 
with very high precision. The numerical solutions are 
further checked using Virial theorem 0. 

To characterize the properties of the equilibrium den- 
sity, we define the aspect ratio of the density profile as 
Ko = \/ {Rx) / {R-l) ■ Since ko = A for a noninteracting 
Fermi gas, kq/A measures the density deformation in- 
duced by dipole-dipole interaction. Figure [1] (a) shows 
the dipolar interaction dependence of kq/X. We imme- 
diately see that the density profile always gets stretched 
along the attractive direction of the dipolar interaction, 
i.e., the axial (radial) direction for positive (negative) 
D. In general, this result does not hold true for a dipo- 
lar condensate, except that the condensate falls into the 
Thomas-Fermi regime due to the large repulsive con- 
tact interaction |l9J. As the dipolar interaction is al- 
ways partially attractive, the system becomes unstable if 
the dipolar interaction exceeds certain critical strength. 
In other words, a dipolar Fermi gas can only be stable 
when D*_ < D < D*^_. To understand how the system 
collapses, we introduce the reduced 'volume' of the cloud 



t^o = 



A-i/2(3/32)i/2' 



which becomes unit in the absence of 
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FIG. 3: (a) The time dependence of cloud aspect ratio for A = 
0.1. In descending order, the dipolar interaction parameters 
are D = —2.3, —1, 0, 1, and 2.4. (b)The dipolar interaction 
dependence of asymptotical aspect ratios for various A's. 



the dipolar interaction. In Fig. [l](b), we plot the D de- 
pendence of i/Q for various trap geometries. Generally, i/q 
increases (decreases) with |-D| if the overall dipolar inter- 
action is repulsive (attractive). For near spherical trap, 
Uq always decreases with |D| when the dipolar interaction 
approaches the critical values. However for highly elon- 
gate or oblate traps, the system becomes unstable even 
when the overall interaction is repulsive. This observa- 
tion suggests that in these cases the collapse is initiated 
locally, which is also confirmed in our numerical simula- 
tion. 

In Fig. [21 we present the A dependence of the criti- 
cal dipolar interaction strength. |1?1| monotonically de- 
creases as one increases A, indicating that a more pro- 
nounced cigar-shaped trap stabilizes the dipolar Fermi 
gas with negative D. This can be easily understood 
as the dipolar interaction is repulsive along z-axis. On 
the other hand, D*^ barely depends on the trap geom- 
etry for A < 20, and it only shows significant increase 
as we further increase A. Compared to the variational 
results Q, the stability gained by increasing A is very 
limited and D'!^ always remains finite. The large discrep- 
ancy between the variational and the numerical predica- 
tions on stability boundaries stems from the fact that a 
simple variational ansatz is generally incapable of captur- 
ing the local collapse. Similar discrepancy between the 
variational and numerical results also exists for dipolar 
Bose-Einstein condensates [1^ [20, [21| . 

Based on above discussion, stabilizing a dipolar Fermi 
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FIG. 4: The D dependence of the release energy for various 
A's. 



gas by increasing A is inefficient. Consider, for exam- 
ple, a typical polar molecule with d = 1 Debye and m = 
lOOa.m.u, and the radial trap frequency uj± = (27r)100 
Hz, one reaches D « 1.5r]N^^^. Therefore, without tun- 
ing 77, a highly oblate trap with A = 100 can only sustain 
about N ~ 4.7 x lO'* fermions. Taking into account the 
deformation of momentum distribution may further lower 
the critical number of molecules [11]. More efficiently, we 
may tune |r/| to a small value to stabilize the system. Al- 
ternatively, we may also load the polar molecules into 
a ID optical lattice, which not only provides an highly 
oblate trap, but also divides the system into many sub- 
systems containing fewer molecules. 

Free expansion. — To study the free expansion, we nu- 
merically evolve Eqs. ([4]) with the restoring force term, 
A|6j, being removed. In Fig. [ni(a), we present the time 
dependence of the cloud aspect ratio Kt = Kobx{t)/bz{t). 
For a noninteracting Fermi gas, Kt always reaches unit 
at the end of expansion, independent of the trap geom- 
etry. In the presence of dipolar interaction, Kt quickly 
reaches its asymptotical value Kqo which deviates from 1 . 
Figure [3l (b) shows the D dependence of the asymptoti- 
cal aspect ratio. One immediately sees that independent 
of the trap geometry and the sign of £), the expanded 
cloud always gets stretched along the attractive direc- 
tion of dipolar interaction. At first sight, this result may 
look counterintuitive, but it can be understood as the 
dipole-dipole mean-field potential accelerates the parti- 
cles along attractive dipolar interaction direction, while 
decelerates those along repulsive direction [l3|. It is also 
worth mentioning that, for polar molecules, the deforma- 
tion in the expanded density should be readily observed 
in experiment. 

In addition to the expanded density profile, the release 
energy i?roh sum of the kinetic and interaction energies 
after switching off the trap, can also be measured exper- 
imentally. The release energy in the absence of dipolar 
interaction is E^J — |(6A)^/'^ for a Fermi gas. The typi- 
cal behavior of the release energy is presented in Fig. [D 
For a highly oblate trap, the release energy may differ 
significantly from that of a noninteracting Fermi gas. 

Deformation of momentum distribution. — As we 




FIG. 5: Plot of contour lines on \-D parameter plane for 
Ecx/Edd = 0.05 (solid line) and 0.1 (dashed line). The direct 
and exchange dipolar interactions ratio for (A, D) enclosed 
by the contour line is lower than that on the corresponding 
contour line. 



have mentioned, the results presented above are obtained 
based on the assumption that the momentum distribu- 
tion is isotropic. Consequently, the exchange dipole- 
dipole interaction vanishes in our treatment. In Ref. 
[ll| . Miyakawa et al. however showed that, due to the 
anisotropy of dipolar interaction, the exchange interac- 
tion may induce the deformation in momentum space, 
such that the momentum distribution is always stretched 
along the attractive direction of the dipolar interaction. 
It is therefore important for us to compare the direct 
and exchange dipolar interactions in order to validate 
our calculations. In Fig. \El we presented the direct 
and exchange dipolar interaction ratio E^d/Eex based 
on the variational calculation developed in Ref. [ll|. It 
can be seen that the exchange interaction has most dra- 
matic effect when the trapping potential is near spheri- 



cal, and it becomes less important in highly prolate and 
oblate traps. The reason behind this is that, compared 
to the density distribution, the momentum distribution 
only weakly depends on the trap geometry. 

The deformation in momentum distribution also gives 
rise to an interesting question about how to observe it 
directly. Since interaction always contaminates momen- 
tum distribution during expansion, we have to turn it off 
during expansion. This can actually be achieved by tun- 
ing 7] to zero simultaneously with the switch-off of the 
external trapping potential, such that the gas expands 
ballistically. 

Conclusion. — We have studied the equilibrium state 
properties and the free expansion of a dipolar Fermi gas 
based on the semi-classical theory. From the stability 
diagram, we propose to stabilize a gas of fermionic po- 
lar molecules by directly tuning the dipolar interaction 
via a fast rotating field. We have also shown that the 
expanded gas always get stretched along the attractive 
direction of the dipolar interaction. For a typical polar 
molecular formed by alkali atoms, the deformation in ex- 
panded density profile can be readily detected. Finally, 
we point out that the deformation of momentum distri- 
bution can be directly observed by switching off the dipo- 
lar interaction during free expansion. We hope that our 
work will stimulate experimental efforts along this line. 
Our future work will include the study of the collective 
excitations of trapped dipolar Fermi gases. 
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